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ABSTRACT 

We present spatial analysis of the deep (57ks) ROSAT HRI X-ray image of the E4 
galaxy NGC 720. The orientation of the HRI surface brightness is consistent with the 
optical position angle (PA) interior to semi-major axis a ~ 60" (optical R e ~ 50"). 
For larger a the isophotes twist and eventually (a ^ 100") orient along a direction 
consistent with the PA measured with the PSPC data (Buote & Canizares 1994) - 
the ~ 30° twist is significant at an estimated 99% confidence level. We argue that this 
twist is not the result of projected foreground and background sources, ram pressure 
effects, or tidal distortions. If spheroidal symmetry and a nearly isothermal hot gas are 
assumed, then the azimuthally averaged radial profile displays features which, when 
combined with the observed PA twist, are inconsistent with the simple assumptions 
that the X-ray emission is due either entirely to hot gas or to the combined emission 
from hot gas and discrete sources. We discuss possible origins of the PA twist and 
radial profile features (e.g., triaxiality) . 



Introduction 



The intrinsic shapes of elliptical galaxies are currently not well understood (e.g., de Zeeuw & 
Franx 1991; de Zeeuw 1995). Hydrodynamic simulations of galaxy halos in a Cold Dark Matter 
universe typically produce flattened (e ~ 0.5), oblate-triaxial halos (e.g., Katz & Gunn 1991; 
Dubinski 1994). Yet there is a paucity of observational constraints for the intrinsic shapes and, 
in particular, the degree of triaxiality in ellipticals. For a handful of galaxies the intrinsic shapes 
have been precisely measured and suggest highly elongated halos with e ~ 0.6 (e.g., Sackett 1995), 
though the degree of triaxiality has not been constrained. Recently, Statler (1994a,b) has shown 
that in principle the stellar velocity fields can effectively probe triaxiality of ellipticals, however the 
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observations required are extensive and the current constraints allow for a large range of triaxial 
shapes. 

The E4 galaxy NGC 720 is a promising candidate for analysis of the shapes of its X-ray 
isophotes because it should not be substantially affected by ram pressure or tidal distortions 
due to it being very isolated from other large galaxies (Schechter 1987). Moreover, NGC 720 is 
quite bright in the ROSAT band (8 x 10 -13 erg cm -2 s _1 for 0.5-2.0 keV) and is likely to be 
dominated by emission from hot gas as indicated by its large ratio of X-ray to optical luminosity 
(Kim, Fabbiano, & Trinchieri 1992). Also, radio emission has not been detected from NGC 720 
(e.g., Birkinshaw & Davies 1985) indicating that the gas should not be substantially disturbed by 
magnetic fields. Since environmental effects should be unimportant, the X-ray-emissivity should 
trace the shape of the three-dimensional gravitational potential, independent of the temperature 
profile of the gas (Buote & Canizares 1994, 1996); though possible rotation of the hot gas may 
affect the emissivity shape in the core (see §|). Thus, the X -ray isophote shapes and orientations 
probe the intrinsic shape of the potential, and hence the mass, of the galaxy. 

From analysis of ROSAT Position Sensitive Proportional Counter (PSPC) data of the E4 
galaxy NGC 720 Buote & Canizares (1994, hereafter BC) concluded that the isophote shapes 
and radial profile of the X-ray emission requires an extended, massive, flattened dark halo 
(e ~ 0.5 — 0.7). However, the X-ray isophotes at a ~ 100" are misaligned by ~ 30° from the optical 
isophotes, a fact that was not incorporated into the analysis of BC. The orientations of the inner 
X-ray isophotes (a <^ 60") are not well determined by the PSPC data because of the width of the 
point spread function (PSF) (~ 30" FWHM). To better understand the nature of this measured 
misalignment between X-ray and optical isophotes, we obtained a high resolution observation to 
specifically address the orientation of the X-ray contours at small radii. 

2. Observations and Data Analysis 

From January 8-12, 1994 NGC 720 was observed with the High Resolution Imager (HRI) on 
board ROSAT (Triimper 1983) resulting in an effective exposure of 57ks; for descriptions of the 
ROSAT X-ray telescope see Aschenbach (1988) and the HRI see David et al. (1995). The spatial 
resolution of the HRI is ~ 4" FWHM, but the HRI possesses only limited spectral resolution in its 
0.1 - 2.4 keV pass-band. 

To prepare the HRI image for spatial analysis we (1) rebinned the image into a more 
manageable 1638 x 1638 field of 2.5" pixels, (2) searched for time intervals where the background 
was anomalously high, (3) flattened the image, (4) selected PHA channels 1-5 to optimize 
signal-to- noise ratio (S/N), (5) examined the accuracy of the aspect solution, (6) removed 
foreground and background sources embedded in the galactic continuum, and (7) subtracted the 
background. All of the reduction procedures were implemented with the standard IRAF-PROS 
and FTOOLS software. 
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By masking out detected sources in the HRI field (including NGC 720 itself using a circle of 
10' radius), we constructed the background light curve which shows no significant fluctuations 
above those expected for the HRI due to external radiation (David et al. 1995). The image was 
flattened with the exposure map generated by the task hriexpmap in FTOOLS; note that the true 
resolution of this exposure map actually corresponds to 5" pixels. The ground calibration results 
of David et al. (1995) show that several PHA channels have low S/N. We examined the effects 
of the PHA channels on S/N for a circle of 50" radius centered on NGC 720; the background was 
taken from an annulus centered on NGC 720 with (r ~ 7', Ar = 1') with sources masked out. 
Including only PHA channels 1-5 optimized the S/N which is consistent with the results of David 
et al., with the possible exception of channel 1 for which David et al. concluded also has low S/N. 

Since we are concerned in this paper with measurements of position angles and ellipticities of 
the X-ray isophotes in NGC 720 an accurate aspect solution for the HRI observation is essential. 
The point sources in the field are not sufficiently bright to usefully perform the aspect correction 
algorithm proposed by Morse (1994). Hence, we examined the aspect solution individually for the 
different observational intervals (OBIs); i.e. those time intervals where the spacecraft continuously 
pointed on NGC 720. From analysis of the point sources in the fields of the three OBIs with 
appreciable exposures we find no statistically significant relative shifts of the OBIs. Typically, the 
required shifts were 0.4 ± 0.8 pixels (i.e. 1" ± 2"). The effects of the small shifts consistent with 
this result are insignificant when considering the measurement of ellipticities and position angles 
on scales of interest to this study (i.e. ^ 15"). 

There is one point source 150" to the SE, another 200" to the NW, and 4 more straddling the 
circle of 400" radius from the galaxy; within ~ 2' of the center of NGC 720 there are no obvious 
point sources. Since we cannot obtain useful constraints on the ellipticity and position angle of the 
X-ray isophotes for r > 100" due to inadequate S/N we simply remove these sources. We removed 
sources by first choosing an annulus around each source to estimate the local background. We 
then fitted a second order polynomial surface to the background and replaced the source with 
the background. This procedure is well suited for estimating quadrupole moments of high S/N 
cluster images; see Buote & Tsai (1996) for a thorough discussion. The effects of unresolved points 
sources on the shapes and orientations of the isophotes are examined in the next section. 

The reduced image is plotted in Figure [j], where for display we have rebinned the image into 
5" pixels and smoothed with a Gaussian of a = 1 pixel. 

The final step in the image reduction is to subtract the background. We only require 
background subtraction for computation of the radial profile, not for computation of the 
ellipticities and position angles. Since we extend the radial profile into regions where the galaxy 
flux <^ background level we need to carefully estimate the background to reduce large fractional 
systematic errors in the radial profile at large radii. The HRI surface brightness profile of NGC 
720 reaches the background level for r ~ 300", somewhat smaller than the outer radius r ~ 400" 
we detected with the more sensitive PSPC (BC). 
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To estimate the background level we compared circular annuli of different radii centered on 
NGC 720. We restricted analysis to regions outside of r = 400" to minimize residual contribution 
from the galaxy. Moreover, to reduce any errors in the HRI vignetting correction and possible 
scattered light in the outer regions of the HRI image we confined our analysis to regions as near 
the center of the field as possible. From consideration of these issues we adopted the 500" — 600" 
annulus to estimate the background level where the residual contribution to the emission from 
NGC 720 should be < 1%. We estimate that this background level is accurate to within ~ 3%. 

3. X-ray Isophote Shapes and Orientations 

As is typical for X-ray images of early-type galaxies the small number of counts (~ 1000) 
for the HRI image of NGC 720 implies that we can only hope to measure with any precision the 
ellipticity and position angle of the aggregate X-ray surface brightness in a large aperture. The 
method we employ is an iterative procedure and is analogous to computing the two-dimensional 
moments of inertia within an elliptical region where the ellipticity, 6m, is given by the square root 
of the ratio of principal moments and the orientation of the principal moments gives the position 
angle, 9m (see BC; Buote & Canizares 1996). The parameters 6m and 6 m are good estimates of 
the ellipticity (e) and position angle (6, or PA) of an intrinsic elliptical distribution of constant 
shape and orientation. For a more complex distribution €m and 6m are average values weighted 
heavily by the outer parts of the regions. 

To estimate the uncertainties on cm and 6m we employ a Monte Carlo procedure described 
by Buote & Tsai (1996) slightly modified for application to HRI images. Because the HRI surface 
brightness profile of NGC 720 appears to be rather complex in shape and orientation, we apply 
the Monte Carlo procedure directly to the HRI image rather than some best-guess model (e.g., 
an elliptical (3 model) which may not account for all of the relevant features. To construct an 
"average" model for the real HRI image we smooth the reduced image (including background) 
with a Gaussian of width a = 2.5" (approximately half the width of the on-axis HRI PSF) to 
better approximate pixels with counts in the real image. We examined the sensitivity of the 
derived values of €m and 6m to the value of a. The estimated 68% confidence limits on em and 
6m for semi-major axes ^> 30" typically vary by < 10% for values of a = 0" — 2.5". For a = 5" 
the derived confidence estimates for these semi-major axes vary by as much as 25%. However, for 
semi-major axes ^ 75" the 68% confidence limits vary by < 5% for a = 0" — 5". 

To the smoothed image we added point sources having spatial properties consistent with the 
HRI PSF and numbers consistent with the logA r (> S) - logS distribution given by Hasinger 
et al. (1993). We excluded from the simulations any bright point sources (i.e. > 5cr above the 
background) appearing within ~ 100" of the galaxy center as they would be easily detected by 
visual examination and removed. Poisson noise was then added to the image consistent with the 
exposure of the NGC 720 HRI observation. 
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We performed 1000 realizations and denned the 90% (for example) confidence limits on cm 
and 6m for each semi-major axis to be the 50th smallest and the 50th largest values obtained from 
the 1000 simulations. Although this definition is arbitrary it provides a simple, realistic measure 
of the significance of the ellipticity and PA. 

In Table [l] we list the results for e m and 9 m and their associated 68% and 90% confidence 
estimates for semi-major axes (a) ranging from 10" to 120". The isophotes are significantly 
elongated with cm ~ 0.3 ± 0.1 for 10" < a < 120". These values for cm are consistent with the 
PSPC values (BC) for a p> 75" but are systematically larger than those obtained with the PSPC 
for a <^ 75", which is consistent with the PSPC values at small a being reduced as a result of its 
larger PSF. 

The position angles of the isophotes vary significantly with semi-major axis; the PA profile 
and the estimated 2a errors are plotted in Figure ^. For a <^ 60", the PA is consistent with the 
optical PA of 142° (e.g., Peletier et al. 1990). At larger a, however, the PA decreases and levels 
off at values consistent with PA = 114°, the position angle obtained from analysis of the PSPC 
at these a. The variation is highly significant: for a = 90" we estimate that PA = 99° — 132° at 
the 95% confidence level, and PA = 90° — 137° at the 99% level, still significantly less than the 
measured value of 144° at a = 60". 

Possible origins for this PA twist are chance alignments of a few appropriately placed 
unresolved foreground and background sources, the residual effects of ram pressure distortions, or 
tidal interactions with neighboring galaxies. For the case of tidal fields, there are no giant galaxies 
within 1 Mpc [Dressier, Schechter, & Rose 1986] which means that the expected gravitational field 
due to a neighbor comparable in size to NGC 720 would contribute only 1% to the potential at 
a radius of 10 kpc from NGC 720. If the other effects are important then the surface brightness 
should display substantial departures from elliptical symmetry. To investigate this possibility we 
first examined whether the centroid of the surface brightness shifts with increasing aperture size. 
Over the whole range of a explored the centroid shifts by less than one pixel in both x and y 
directions; i.e. there is no significant centroid shift. 

We also examined higher order symmetry. First, we aligned coordinate axes centered on the 
galaxy with the x axis at the optical PA = 142°. Within a radius of r = 60" we computed the 
counts in the four quadrants in these coordinates. The results are listed in Table The counts 
are consistent with a mean value of ~ 285 and their la Poisson uncertainties of 22; i.e. within 
r = 60" there is no evidence for asymmetry between the four quadrants. By instead aligning 
the coordinate axes with PA = 114° and examining the surface brightness in an annulus from 
r = 60" — 90" we explored the symmetry in the outer regions; see Table |. A gain, there is no 
evidence for asymmetries. 

Unresolved point sources (and statistical noise) cannot produce the observed PA twist as 
demonstrated by the above Monte Carlo simulations (Table |l|). To further investigate this issue 
we estimated the probability that unresolved sources contribute sufficiently to (specifically) the 
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outer, low surface brightness, isophotes to produce the position angle twist. The region of interest 
lies between the isophotes with a ~ 75" — 100". If part of this region is significantly contaminated 
by point-source emission, then a PA twist could result. Consider, e.g., a 25" x 25" square within 
this region in which the flux from NGC 720 is measured to be ~ 4 x 10 -15 erg cm -2 s _1 in the 
0.1-2.4 keV band using the PSPC image. From the log N(> S) - log S distribution for ROSAT 
X-ray sources (Hasinger et al. 1993) we estimate the probability that there is a source with flux 
at least half of the measured value in the square region to be only ~ 2%. Hence, it is unlikely 
that the gross features of the surface brightness in the key regions r ~ 75" — 100" are affected 
substantially by the emission from unresolved foreground and background sources. 

4. Radial Profile 

Using the background-subtracted HRI image we constructed the azimuthally averaged radial 
profile located at the centroid of the galaxy emission. The centroid was computed in a circular 
aperture of radius 120" which effectively encloses ~ 85% of the total flux of the galaxy; note 
that the centroid varies by < 2.5" for smaller aperture sizes. The radial profile was rebinned 
so that each bin had S/N > 6. This yields six 5" bins from r = 0" — 30", three 15" bins from 
r = 30" — 75", and one 45" bin from r = 75" — 120". The radial profile is plotted in Figure |3[ 

For r ;> 30", the radial profile is consistent with a power-law form. However, the profile 
appears to level off between r ~ 10" — 30" and then rises again for r <^ 10". These qualitative 
features of the radial profile appear to be robust to small changes in the sizes of the radial bins. 
As we show in the next section, although these features are not of sufficient S/N to distinguish 
between a wide range of models, they are of sufficient quality to eliminate some of the simplist 
models that previously were consistent with the X-ray data for NGC 720. 

5. Analysis and Discussion 

Since the emission from unresolved foreground and background sources and effects from ram 
pressure and tidal distortions appear to be insignificant (see §||), we now consider origins of the 
surface brightness features that are intrinsic to the galaxy. First we examine whether the HRI data 
may be described by a model where the X-ray emission is entirely due to hot gas in hydrostatic 
equilibrium with the underlying potential of NGC 720. We shall consider two classes of models - 
(3 models and power-law or Hernquist Spheroidal Mass Distributions (SMDs). First we take the 
surface brightness distribution to be given by the (3 model (Cavaliere & Fusco-Femiano 1976), 
Yih g (R) = So[l + (i?/-R c ) 2 ]~ 3 ^ +1 / 2 , which is a power law for R ^> R c . This model corresponds to a 
logarithmic potential at large radii (e.g., Appendix C. of Trinchieri et al. 1986) and is generally 
a good fit to the X-ray radial profiles of ellipticals observed with the Einstein IPC (e.g., Forman, 
Jones, & Tucker 1985; Trinchieri et al. 1986) and the ROSAT PSPC - in particular that of NGC 
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720 (see BC). 

To perform a more consistent comparison between the model and data we evaluate the model 
on a grid of pixels corresponding to the same scale as the data and convolve the model with the 
relevant PSF. The radial profile is then constructed and binned in the same manner as done for 
the data. This procedure does yield slightly different parameter values and x 2 values than if the 
convolved model is directly compared to the radial profile of the data as is typically done. In all 
of the fits we set the ellipticity of the hot gas to e x = 0.25 which is essentially the mean value 
of the X-ray isophotes. The fitted parameters and \ 2 values change by £ 5% when considering 
ellipticities e x = 0.2 — 0.3. 

The (3 model is a poor fit to the HRI radial profile: x 2 = 18 for 7 dof. The poor quality of 
the fit is due entirely to the inner radial bins (r < 30"). In contrast, BC showed that the [3 model 
provides a good fit to the PSPC profile of NGC 720. We reanalyzed the PSPC data and rebinned 
the PSPC radial profile so that each bin has S/N > 5 and energy range 0.5-2.0 keV; BC instead 
fixed the bin size to 15". The (3 model applied to this rebinned PSPC radial profile yields a 
formally marginal fit: x 2 = 19 for 10 dof. Similar to the HRI data, the major source of discrepancy 
between model and data is due to the two inner bins (0" — 15", 15" — 30") - the bin near 90" is not 
as important due to its larger fractional statistical uncertainty. The best-fit parameters for the 
HRI and PSPC fits are listed in Table ||. The HRI profile is shallower (i.e. smaller ff) than that 
of the PSPC which is probably an artifact of the weak constraints at large radii (r p> 100"); i.e. 
when the PSPC data is fit only to bins within r = 100" then smaller (3 values are also obtained. 

It is necessary to fit the model jointly to the HRI and PSPC data in order to obtain constraints 
acceptable for both data sets; the model for each data set is convolved with its appropriate PSF. 
By doing this we now have four free parameters: a x , (3, and the normalizations for the HRI and 
PSPC radial profiles. The joint fit is shown in Figure [3| and the best-fit and 90% confidence levels 
(2 interesting parameters) are shown in Table |3|. The fit is formally worse ( x 2 = 53 for 19 dof) 
than the fits to the data sets individually which is primarily due to joining the shallower HRI 
profile to the steeper PSPC profile; note the poor fit is not the result of the small (< few %) 
uncertainties in the background estimates. 

These qualitative and quantitative results are reproduced when we consider the spheroidal 
mass models (SMDs) used in BC that are consistent with the PSPC spatial and spectral (i.e. 
nearly isothermal) data; i.e. models where the mass density is stratified on concentric, similar 
spheroids and has either a power-law or Hernquist (1990) radial profile. For example, a model 
with mass density p ~ r~ 2 and an ellipticity of 0.50 fit jointly to the HRI and PSPC data yields 
a poor fit: \ 2 = 46 for 19 dof. Similar results are obtained for Hernquist (1990) models and its 
generalization to shallower profiles (i.e. p ~ r _1 (a c + r) _1 ). Hence, simple plausible models for the 
hot gas cannot fit the radial profiles of the HRI and PSPC data, where the disagreement between 
models and data primarily occurs in the inner regions r <J 30". 

We are thus compelled to increase the complexity of our model to fit the X-ray surface 
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brightness. There must be significant contribution to the X-ray emission from discrete sources 
in NGC 720; i.e., the integrated emission from X-ray binaries (see, e.g., Canizares, Fabbiano, & 
Trinchieri 1987; Kim, Fabbiano, & Trinchieri 1992). This "discrete" component should have a 
spectrum that is harder than that of the hot gas and essentially similar to that found in spirals 
where discrete sources are likely to account for most of the emission. 

Considering the relatively large value of Lx/Lb for NGC 720 we do not expect a large 
contribution to the X-ray emission from discrete sources in the 0.1 - 2.4 keV energy range (Kim et 
al. 1992). We estimate the relative flux of the hot gas and a harder component using the PSPC 
spectrum computed in a circle of 120" about the galaxy centroid. We fitted the spectrum with 
two Raymond-Smith (1977, updated to current version) models where we fixed the column density 
to its Galactic value (Stark et al. 1992), the abundances to 20% solar (approximately the best-fit 
single-component value - see BC), and the temperature of the hard component to 4.5 keV which 
corresponds to the integrated spectrum of Sb+Sc galaxies analyzed by Kim et al. (1992). This 
temperature is also consistent with the hard component determined for NGC 4472 with ASCA 
(Awaki et al. 1994; Matsushita et al. 1994), a galaxy that has a value of Lx/Lb very similar to 
NGC 720. We find for these models that the emission of the hard component can make up at 
most 32% of the total emission within 120" (90% confidence on two interesting parameters) but 
with a lower limit approaching zero (0.1%). This range is consistent with the ASCA observation 
of NGC 4472 where Awaki et al. (1994) determined that only ~ 10% of the X-ray emission in the 
0.5 - 5 keV band is due to a hard component (for about twice the radius used for NGC 720). In 
the 0.1-2.4 keV band of the HRI and PSPC this percentage should be even smaller. 

We place alternative constraints on the total emission from discrete sources using the 
spatial properties of the X-ray surface brightness. Although there exists no strong observational 
constraints on the spatial distribution of such a component in ellipticals, we make the reasonable 
assumption that the discrete emission has a spatial distribution similar to that of the optical 
light. We take the surface brightness of the discrete component, E^, to have the two-dimensional 
distribution of the -B-band light of NGC 720: the HST data from Lauer et al. (1995) define the 
core while the rest of the galaxy is defined by the data from Peletier et al. (1990), who find no 
evidence for color gradients in NGC 720. The hot-gas component, E^ s , is again modeled as either 
the (5 or the isothermal power-law or Hernquist (1990) SMDs. 

In order to specify the parameters of this discrete + hot gas model we fitted this model jointly 
to radial profiles of the HRI and PSPC data, convolved with the appropriate PSFs, and using the 
same procedure as above. The HRI data provide the dominant constraint at small radii (R <^ 30") 
while the PSPC data dominates over the rest of the galaxy (R ~ 30" — 400"). For each fit the 
ratio of emission from hot gas to the emission from discrete sources, E^/E^, is computed in a 
circle of radius 120" from the galaxy centroid. 

The results of the joint fits to the HRI and PSPC radial profiles are listed in Table Q and 
shown in Figure [3|. Let us first consider the case where the hot gas is taken to be a (5 model. The 
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range of ratios of discrete emission to hot gas obtained are consistent with the above spectral 
results from the PSPC which serves as a comforting consistency check between the spatial and 
spectral methods. Although only one free parameter has been added (S^/S^) to our model, 
the quality of the fit is substantially improved (x 2 = 29 for 18 dof), though still formally only 
a marginal fit (P x 2 = 95%). Again the discrepancy between model and data may be attributed 
primarily to the inner 30" region. The model fits the PSPC profile acceptably, but the HRI profile 
(particularly the core) is poorly fit: Xpspc = ^ f° r ^ dof and Xhri = 15 for 6 dof. In contrast to 
the hot gas alone, adding the discrete component allows the innermost bin to be well fitted (i.e. 
this combined hot gas + discrete model can account well for the excess emission in the central 
bin). However, this model does not produce the flattening of the HRI radial profile between 
r ~ 10" — 30", which is the dominant source of discrepancy between model and data. When the 
hot gas is modeled using the SMDs (see above) we obtain results in excellent agreement with the 
(5 models (see Table ||) . 

Although the formal fit to the data is marginal, it is worth examining whether the discrete 
+ hot gas model can reproduce the observed position angle twist. We take the discrete model 
to be oriented along the optical PA while the hot gas is oriented along the PA of the outer 
isophotes {PA = 114°). We find that the models do not produce large PA twists, the largest 
twist is generally produced by the best-fit model to the radial profile where APA ~ 10° between 
semi-major axes of 60" and 100". The statistical uncertainties on simulated observations of these 
models are rather large and similar to what was found for the data (Table |l]): the observed PA 
twist between semi-major axes of 60" and 100" is inconsistent at the estimated ~ 90% confidence 
level. Since, however, an acceptable model must produce the PA profile between 60" and 100" (not 
simply the endpoints), the discrepancy may be more pronounced. (This issue requires significantly 
more computational effort to properly address than is appropriate for our present study.) Hence, 
the discrete + hot gas model cannot fully produce the observed large PA twist (though it will 
contribute) or the features in the core of the radial profile. 

Joint consideration of the HRI and PSPC data of NGC 720 suggests that the X-ray emission 
is more complex than simple isothermal spheroidal models of the hot gas combined with emission 
from discrete sources. Detailed investigation of more sophisticated models is beyond the scope of 
this paper, but we briefly describe some possible solutions. 

Triaxiality: A PA twist can be produced by the projection of a triaxial distribution whose 
three-dimensional isodensity surfaces have shapes which vary as a function of radius (e.g., Mihalas 
& Binney 1981). The ellipticities of the optical isophotes do not vary substantially over the region 
of interest (a ^ 30"; e.g., Peletier et al. 1990). However, the stars do not trace the same shape 
as the potential (because of anisotropic velocity dispersions) and thus do not in general have the 
same three-dimensional spatial distribution, and thus projected distribution, as the gas. N-body 
simulations of galaxies and clusters generically produce strongly triaxial halos (e.g., Frenk et al. 
1988; Dubinski &: Carlberg 1991). In fact, hydrodynamic simulations often produce halos that 
are nearly oblate in the central regions while the outer parts, particularly for clusters formed 
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in filaments, are nearly prolate (e.g., Katz & Gunn 1992; Dubinski 1994; Buote & Tsai 1995). 
Perhaps a large observed PA twist can result if the transition from oblate (r <J 60") to prolate 
(r ;> 100") is sufficiently rapid. However, using simple models with the above specifications and 
a power-law variation in axial ratio between a ~ 60" and a ~ 100" we find it difficult to produce 
PA twists larger than ~ 15 — 20° which also produce the observed ellipticities of the isophotes. 

Intrinsic Misalignment of Light and DM Halo: If the principal axes of the stars and dark 
matter halo are intrinsically misaligned then a PA twist would naturally occur. Such a scenario 
has been proposed by Ostriker & Binney (1989) to explain the origin of warps in disk galaxies. 
However, Dubinski & Kuijken (1995) have shown that large misalignments will rapidly decay 
because the disk, with its large spin misaligned with the halo, will precess and suffer dynamical 
friction with the halo, and thus settle into the halo equatorial plane; often the inner halo actually 
aligns with the disk because the disk is difficult to re-orient due to its large spin. Because a 
giant elliptical like NGC 720 does not rotate substantially (v/a* = 0.15, e.g., Fried & Illingworth 
1994), it is not clear whether these arguments imply the same short timescales. Perhaps 
another mechanism (e.g., phase mixing) ensures alignment for an elliptical. Even if the timescale 
for alignment is found to be short compared to a Hubble time (thus ruling out a primordial 
misalignment), it may be that the misalignment is caused by a dwarf companion that has recently 
struck NGC 720; there are several known dwarf companion galaxies around NGC 720 (Dressier et. 
al. 1986). 

Rotating Cooling Flow: The flattening of the radial profile between r ~ 10" — 30", though not 
of high enough significance to warrant detailed modeling, deserves some attention. It is possible 
that the surface brightness features on these small scales are due to patchy absorption by cold 
clouds condensed out of the hot phase (e.g., Fabian 1994). However, another aspect of the cooling 
gas may be the origin for the features. Even though giant ellipticals like NGC 720 typically rotate 
slowly, in principle slowly rotating gas at large galacto-centric radii can cool, flow towards the 
center, eventually be driven to dynamically important rotational velocities if angular momentum 
is conserved (e.g., Cowie, Fabian, & Nulsen 1980), and thereby form a spinning disk (e.g., Kley & 
Mathews 1995). Not only will the ellipticity of the X-ray isophotes become more elongated in the 
region of dynamically important rotation, the shape of the radial profile will be affected as well 
(e.g., Kley & Mathews 1995). 

6. Conclusions 

The interesting features of the HRI surface brightness of NGC 720 provide more questions 
than answers into the origin of its X-ray emission and the distribution of its underlying mass. 
Presently ASCA, with its improved spectral resolution over ROSAT, can help elucidate the role 
of the total emission from discrete sources in NGC 720. However, we must await the launch of 
AXAF to have the combined spatial and spectral resolution to distinguish models of the origin of 
the position-angle twist and flattening of the radial profile of the HRI data. 
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Table 1: X-ray Ellipticities and Position Angles 



a 

(arcsec) 


£M 


68 ( 


Yo 


90% 




68 ( 


Vo 


90% 


Counts 


10 


0.36 


0.06 - 


0.37 


0.02 - 


0.51 


75 


45 - 


132 


19 - 


160 


158 


20 


0.14 


0.08 - 


0.28 


0.03 - 


0.36 


119 


84 - 


144 


60 - 


159 


418 


30 


0.39 


0.22 - 


0.39 


0.16 - 


0.44 


129 


131 - 


156 


123 - 


165 


527 


40 


0.35 


0.25 - 


0.40 


0.20 - 


0.45 


155 


145 - 


162 


139 - 


167 


720 


50 


0.28 


0.25 - 


0.37 


0.20 - 


0.41 


151 


144 - 


157 


139 - 


162 


908 


60 


0.35 


0.23 - 


0.38 


0.17 - 


0.40 


144 


139 - 


153 


133 - 


158 


977 


70 


0.30 


0.24 - 


0.38 


0.18 - 


0.44 


134 


126 - 


144 


120 - 


149 


1130 


80 


0.33 


0.25 - 


0.39 


0.20 - 


0.42 


124 


114 - 


130 


110 - 


135 


1220 


90 


0.28 


0.23 - 


0.36 


0.18 - 


0.41 


118 


108 - 


125 


102 - 


130 


1351 


100 


0.28 


0.17 - 


0.34 


0.12 - 


0.38 


112 


103 - 


123 


95 - 


129 


1427 


110 


0.24 


0.19 - 


0.35 


0.12 - 


0.39 


110 


99 - 


124 


87 - 


131 


1537 


120 


0.28 


0.18 - 


0.35 


0.09 - 


0.40 


111 


97- 


121 


89 - 


130 


1611 



Note. — The values of em (and confidence limits) are computed within an aperture of semi-major axis a on the 
image having 2.5" pixels with the background included; the counts, however, have the background subtracted. 



Table 2: Symmetry of X-ray Surface Brightness 





PA = 142° 


PA = 114° 


Quadrant 


(r = 60") 


(60" < r < 90") 


I 


303 ± 23 


86 ± 18 


II 


253 ± 21 


96 ± 18 


III 


287 ± 22 


82 ± 18 


IV 


284 ± 22 


99 ±18 



Note. — Background-subtracted counts and their associated ler uncertainties assuming Poisson statistics. See 
for an explanation of the regions. 
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Table 3: Radial Profile Models 





R c 










Model 


(arcsec) 




^ds/^hg 


x 2 


dof 


(3- HRI 


6±3 


0.40 ± 0.03 




18 


7 


p - PSPC 


18 ±4 


0.49 ± 0.02 




19 


10 


(3 - BOTH 


13 ±3 


0.46 ± 0.02 




53 


19 


(3 + DISCRETE 


3O+1 


o.5ol8:8f 


o.4ot8:I5 


29 


18 


SMD - BOTH 








46 


19 


SMD + DISCRETE 






o.4o+8:I? 


29 


18 



Note. — Under "Models" , /3 stands for the hot gas represented by the (3 model while SMD indicates the spheroidal 
mass models found by BC to be consistent with the PSPC data for NGC 720. "BOTH" indicates model is fit jointly 
to HRI and PSPC data. The DISCRETE models are fit jointly to both data sets. For the SMD models we only list 
the value of \ 2 and the flux ratio for the DISCRETE case. 
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Fig. 1.- 

Contour map of the HRI X-ray surface brightness of the E4 galaxy NGC 720 binned into 5" pixels 
for display; the contours are separated by a factor of 2 in intensity and the direction of Celestial 
North is up and East to the left. The image has been corrected for the effects of exposure 
variations and telescopic vignetting. The image has been smoothed for visual clarity with a 
Gaussian of a = 1 pixel, although the image used for analysis is not smoothed as such. 




NGC 720 



Fig. 2.- 

The HRI position angle computed within a specified semi-major axis, a, and the estimated 95% 
confidence intervals; note the position angles are correlated since each a contains the regions < a. 
The optical and PSPC PA are shown approximately over their relevant range of validity. 
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Fig. 3.- 

(a) The HRI background-subtracted radial profile binned such that each bin has S/N > 6. The 
solid line represents the best-fit (3 model and the dotted line corresponds to the best-fit (3 + 
discrete model. The midpoint of each bin is used to define its location in the plot; note the lines 
are only defined at the midpoints of the bins. The dashed horizontal lines represent the 
background levels, (b) The PSPC background-subtracted radial profile binned such that each bin 
has S/N > 5. The squares and crosses represent the corresponding best-fit models from (a). 
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